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ABSTRACT 


In  a  typical  self-contained  space  navigation  system  celestial  observa¬ 
tion  data  are  gathered  and  processed  to  produce  estimated  velocity  corrections. 
The  results  of  this  paper  provide  a  basis  for  determining  the  best  celestial 
measurements  and  the  proper  times  to  implement  velocity  corrections. 

Fundamental  to  the  navigation  system  is  a  procedure  for  processing 
celestial  measurement  data  which  permits  incorporation  of  each  individual 
measurement  as  it  is  made  to  provide  an  improved  estimate  of  position  and 
velocity.  In  order  to  "optimize"  the  navigation,  a  statistical  evaluation  of  a 
number  of  alternative  courses  of  action  is  made.  The  various  alternatives, 
which  form  the  basis  of  a  decision  process,  concern  the  following: 

1.  Which  star  and  planet  combination  provide  the  "best"  available 
observation  ? 

2.  Does  the  best  observation  give  a  sufficient  reduction  in  the  pre¬ 
dicted  target  error  to  warrant  making  the  measurement? 

3.  Is  the  uncertainty  in  the  indicated  velocity  correction  a  small  en¬ 
ough  percentage  of  the  correction  itself  to  justify  an  engine  re¬ 
start  and  propellant  expenditure? 

Numerical  results  are  presented  which  illustrate  the  effectiveness  of 
this  approach  to  the  space  navigation  problem. 
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1.  INTRODUCTION 


During  the  past  two  years,  the  problems  of  guiding  a  space  vehicle 

during  the  midcourse  phase  of  its  mission  have  been  extensively  explored  at 

the  MIT  Instrumentation  Laboratory.  Following  the  specific  demonstration 

of  the  technical  feasibility  of  an  unmanned  photographic  reconnaissance  flight 

to  the  planet  Mars  reported  by  Laning,  Frey,  and  Trageser^,  the  detailed 

navigational  aspects  of  such  a  venture  were  developed  ^  by  Dr.  J.  H.  Laning, 

Jr.,  and  the  present  author.  Later,  a  variable  time  of  arrival  navigation  theory 
(3) 

was  devised  and  contrasted  with  the  earlier  fixed  time  of  arrival  scheme. 
More  recently,  the  question  of  optimum  utilization  of  navigation  data  has  been 
given  considerable  study.  It  is  the  solution  of  this  problem  which  forms  the 
subject  of  the  present  paper. 

The  general  method  of  navigation  is  based  on  perturbation  theory  so 
that  only  deviations  in  position  and  velocity  from  a  reference  path  are  utilized. 
Data  is  gathered  by  an  optical  angle  measuring  device  and  processed  by  a 
spacecraft' digital  computer.  Periodically,  small  changes  in  the  spacecraft 
velocity  are  implemented  by  a  propulsion  system  as  directed  by  the  computer. 

Basically,  three  problems  are  considered  in  this  paper:  (l)  to  identify 
the  best  sources  of  data  available  to  the  space  vehicle  navigator;  (2)  to  define 
the  optimum  linear  operations  for  processing  the  data  in  a  manner  consistent 
with  the  mission  objectives;  and  (3)  to  minimize  both  the  amount  of  navigational 
data  and  the  number  of  corrective  maneuvers  required  without  unduly  comprom¬ 
ising  mission  accuracy. 

The  formulation  of  an  optimum  linear  estimator  as  a  recursion  oper¬ 
ation  in  which  the  current  best  estimate  is  combined  with  newly  acquired 

( 4 ) 

information  to  produce  a  still  better  estimate  was  presented  by  Kalman 
The  author  is  indebted  to  Dr.  Stanley  F.  Schmidt  for  directing  his  attention  to 
Kalman's  excellent  work.  In  fact,  the  original  application  of  Kalman's  theory 

(5) 

to  space  navigation  was  made  by  Schmidt  1  and  his  associates. 
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The  work  described  in  the  following  sections  of  this  paper  was  done 
without  any  detailed  knowledge  of  Schmidt's  activities.  As  a  result  of  this  in¬ 
dependent  approach,  several  new  and  interesting  ideas  have  developed.  Spec¬ 
ifically,  an  extremely  simple  derivation  of  the  optimum  linear  operator  has 
been  achieved  using  only  the  basic  technique  of  least  squares  estimation.  In 
addition,  the  following  new  results  are  noted:  (1)  the  mean- squared  velocity 
correction  is  expressed  directly  in  terms  of  initial  orbital  injection  errors, 
the  errors  associated  with  navigational  measurements,  and  the  errors  in 
establishing  the  desired  velocity  correction  --  hence,  a  statistical  simulation 
of  the  navigation  scheme  may  be  made  without  resorting  to  Monte  Carlo 
techniques;  (2)  a  detailed  procedure  for  incorporating  cross-correlation 
effects  of  random  measurement  errors  in  determining  the  optimum  linear 
operator  has  been  developed;  and  (3)  the  mathematical  problem  of  determining 
the  optimum  plane  in  which  to  make  a  star-planet  angular  measurement  has 
been  solved. 

Throughout  the  paper,  we  shall  deal  exclusively  with  discrete  infor¬ 
mation;  observations  or  velocity  corrections  are  made  at  specific  points  in 
time  which  are  termed  "decision  points."  The  interval  between  decision 
points  is  not  necessarily  uniform  and  may  be  selected  somewhat  arbitrarily; 
e.g.,  the  interval  length  required  for  accurate  numerical  integration  of  the 
trajectory  equations  was  used  in  preparing  the  computational  data  presented 
in  Section  6. 

Finally,  a  few  remarks  relevant  to  notational  conventions  are  appro¬ 
priate.  We  shall  deal  generally  with  both  three-  and  six-dimensional  vectors. 
A  column  vector  of  any  dimension  is  represented  by  a  lower  case  underscored 
letter.  Matrices  are  denoted  by  capital  letters  and  can  be  either  square  or 
rectangular  arrays.  The  transpose  of  a  vector  or  a  matrix  will  be  denoted  by 

a  superscript  T.  Thus,  the  scalar  product  of  two  vectors  a  and  b  will  be 
T 

written  as  a  b.  In  like  manner  a  quadratic  form  associated  with  a  square 

T 

matrix  A  will  be  written  as  jc  Ax,  Finally,  the  expected  value  of  a  random 
vector  x  will  be  indicated  by  an  overscore;  thus,  x  denotes  the  average 
value  of  x. 
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2.  OUTLINE  OF  THE  NAVIGATION  AND  GUIDANCE  PROCEDURE 


2.1  A  Deterministic  Method 

The  basic  process  involved  in  determining  spacecraft  position  by  means 
of  a  celestial  fix  consists  fundamentally  of  a  sequence  of  measurements  of  the 
angles  between  selected  pairs  of  celestial  objects.  Three  independent  and 
precise  angular  measurements  made  at  a  known  instant  of  time  suffice  to 
determine  uniquely  the  position  of  the  vehicle.  Practical  constraints,  however, 
preclude  simultaneous  measurements  without  severely  complicating  the  instru¬ 
mentation.  On  the  other  hand,  if  the  vehicle  dynamics  are  governed  by  known 
laws  and  if  deviations  from  a  pre-determined  reference  trajectory  are  kept 
sufficiently  small  to  permit  a  linearization  of  the  navigation  problem,  then  the 
question  of  simultaneous  measurements  loses  it  significance. 


Under  the  assumptions  of  a  linearized  theory,  a  single  observation 
serves  to  fix  the  position  of  the  spacecraft  in  one  coordinate.  For  example, 
if  An  is  the  angle  measured  at  time  tn  and  is  defined  by  the  lines-of-sight 
from  the  vehicle  to  a  star  and  to  a  nearby  celestial  body,  the  position  of  the 
vehicle  is  established  along  a  line  normal  to  the  direction  toward  the  near  body 
and  in  the  plane  of  the  measurement.  It  is  shown  in  Appendix  A  that  the  devi¬ 
ation  in  position  6rn  of  the  spacecraft  from  the  reference  position  is  related 
to  the  deviation  in  angular  measurement  6A^  by 


8  A  -h.  Sr 

n  -n  — n 

if  the  observation  is  made  at  a  known  instant  of  time  t 


(2.  1) 


The  vector  h  dep- 
-n  r 


ends  upon  the  geometrical  configuration  of  the  relevant  celestial  objects  at 

time  t  as  well  as  the  type  of  measurement  made. 

n  J  r 


Because  of  the  inherent  dynamic  coupling  of  position  and  velocity,  the 

result  at  a  later  time  t  .  ,  of  a  measurement  made  at  time  t  does  not  lend 

n+1  n 

itself  to  simple  geometric  interpretation.  In  order  to  provide  a  geometrical 
description,  it  is  convenient  to  introduce  the  concept  of  a  six  dimensional 
space  in  which  the  coordinates  represent  the  components  of  both  position  and 
velocity  deviations  of  the  vehicle  from  the  reference  path  as  functions  of  time. 


<* 


L 
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Points  in  this  space  are  defined  by  the  six  dimensional  deviation  vector 


8*n  = 


8-rn 
8  v. 


(2.2) 


where  6yn  is  the  deviation  in  the  vector  velocity  of  the  vehicle  from  the  refer¬ 
ence  value.  The  vector  6xn  defines  the  "state"  of  the  vehicle  dynamics  at 
time  t  .  Transition  from  one  state  to  another  is  provided  by  the  matrix  oper¬ 
ation 


n  +  l,n  ^  1  *n  +  1' 


which  is  frequently  referred  to  as  the  "transition  matrix". 

ship  between  6x  and  6x  is  simply 
-n+1  -n  J 

8-xn  +  l  =  ^n+l.n  8*n 

as  shown  in  Section  3.4. 


Indeed,  the  relation- 
(2.  3) 


By  means  of  the  rectangular  matrix  K  defined  by 

K  -  1 

K  0 

Eq.  (2.  1)  may  be  written  in  terms  of  fix  as 

~n 


8  A  =  hT  KT 

n  —  n 


8  -n 


(2.4) 


(2.  5) 


The  submatrices  I  and  O  are,  respectively,  the  three  dimensional  identity 
and  zero  matrices.  Now,  by  combining  Eqs.  (2.3)  and  (2.  5) 

8An^blKT4>nT+1,n8xn  +  1  (2.6) 

it  is  clear  that  the  effect  at  time  t  .  .  of  an  observation  at  time  t  is  to  determine 

n+1  n 

the  component  of  the  six  dimensional  deviation  vector  in  the  direction  defined  by 
T- 1 

the  vector  <i>  .  ,  Kh  .  Six  observations  made  at  different  times  would 

n+ 1 ,  n  —  n 

provide  a  set  of  six  equations  of  the  form  of  Eq.  (2.  6)  .  If  no  two  of  the  com¬ 
ponent  directions  were  parallel,  then  the  deviation  vector  could  be  obtained  by 
inverting  the  six  dimensional  coefficient  matrix. 


2.  2  Statistical  Parameters  of  the  Navigation  Problem 

Because  of  the  presence  of  instrument  inaccuracies  additional  observations 
may  be  used  to  reduce  the  errors  associated  with  the  simple  deterministic  pro¬ 
cess  just  described.  By  applying  least  square  techniques  to  the  observed  data, 
a  more  accurate  estimate  of  position  and  velocity  is  frequently  possible  than 
could  be  obtained  from  the  minimum  number  of  measurements.  For  this 
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purpose,  it  is  necessary  to  know  certain  statistical  information  with  respect 
to  the  instrument  inaccuracies.  In  a  linear  least  squares  estimation  procedure 
all  statistical  calculations  are  based  on  first  and  second  order  averages  and  no 
additional  statistical  data  is  needed. 


At  this  point  of  the  discussion  it  is  necessary  to  distinguish  measured 

fV 

values,  estimated  values  and  true  values  of  various  quantities;  e.  g. ,  6A^  will 

be  the  measured  value  of  the  deviation  in  the  angle  A^  from  its  reference  value 

at  time  t  ,  6A  the  true  value  of  the  deviation,  and  6  A  the  estimated  value.  If 
n'  n 

we  write  ^ 

SAn=SAn+a„  (2.7) 

then  an  will  be  the  error  in  the  measurement.  In  the  subsequent  analysis 

will  be  regarded  as  a  random  variable  with  an  average  value  and  a  variance 


(2.  8) 


The  possibility  of  cross-correlation  of  measurement  errors  will  not  be  excluded; 


i.e.,  in  general,  the  average  unam  may  be  different  from  zero. 


In  Section  4  an  estimation  procedure  is  developed  for  determining  an 
optimal  linear  estimate  of  denoted  bs  6  x^.  As  each  measurement  is  made, 

the  estimate  6x^  is  updated  by  a  simple  recursive  formula  and,  thereby,  the  prob¬ 
lem  associated  with  inverting  sixth  order  matrices  is  avoided.  An  integral 
part  of  the  estimation  technique  is  the  correlation  matrix  of  the  errors  in  the 
estimate.  If  we  write 

S?n  =  ^ -n  +®n  <2-9) 

then 

(2.  10) 

is  the  six  dimensional  error  vector  and  may  be  partitioned  as  shown  using 
e  n  and  6  to  denote,  respectively,  the  position  and  velocity  errors.  The  corr¬ 
elation  matrix  is  thus  defined  by 


(2.  11) 
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When  cross-correlation  of  measurement  errors  is  considered,  a  corr¬ 
elation  vector 


-?nm  an-m  (2.  12  ) 

is  needed  to  represent  the  average  value  of  the  product  of  a  measurement  error 

at  time  t  and  the  estimation  error  at  time  t  . 

n  m 


It  is  important  to  distinguish  between  a  new  estimate  Sx^,  obtained  by 
incorporating  an  observation  at  time  t  ,  and  an  estimate  simply  extrapolated 
from  a  previous  estimate.  For  the  latter  case,  the  notation  is  used  where 

S*n  =<V  n-  1  S?n-  1  (2.  13) 

In  like  manner,  we  define  an  extrapolated  error  vector  e'  and  a  cross- corr¬ 
elation  vector  6  '  =  a  e'  The  extrapolated  correlation  matrix  is 

—  nm  n  -  m 

readily  shown  to  be 


E'  =  <J> 

n  n,n- 


1  En 


*T  1 
n,  n  -  1 


(2.  14) 


As  a  final  comment,  note  that  an  estimate  of  the  deviation  in  the  angle  to  be 
measured  at  time  t  may  be  obtained  from  the  extrapolated  estimate  of  6x^ 
We  have 

=  hj  KT  Sx'  (2.  15) 

and  it  is  this  quantity,  compared  with  the  measured  deviation  6  A^,  which  is 
used  in  arriving  at  a  revised  estimate  of  6xn- 

2.  3  Summary  of  the  Navigation  and  Guidance  Equations 


In  the  navigation  and  guidance  theory  presented  here,  the  problem  of 
launch  guidance  from  Earth  is  ignored.  It  is  assumed  that  the  main  propulsion 
stages  are  completed  at  time  t^  and  that  the  correlation  matrix  Eq  =  E(t^)  is 
specified  initially  from  a  statistical  knowledge  of  injection  guidance  errors. 
The  initial  estimate  of  position  and  velocity  deviation  6x  =  6  x  (tT  )  is  zero, 
since,  in  the  absence  of  any  observation,  the  best  unbiased  estimate  is  that  the 
spacecraft  is  on  course. 


The  time  interval  from  launch  to  arrival  time  t^  at  the  target  point  is 
considered  to  be  subdivided  into  a  number  of  smaller  intervals  by  the  sequence 
of  times  t^,  t^,  .  .  .  called  "decision  points".  At  each  decision  point  one  of 
three  possible  courses  of  action  is  followed:  (1)  a  single  observation  is  made; 
(2)  a  velocity  correction  is  implemented;  or  (3)  no  action  is  taken.  A  revised 
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estimate  of  the  deviation  vector  6x(t)  is  made  at  each  such  point  --  the  form  of 
the  revision  depending,  of  course,  on  the  nature  of  the  decision.  Specifically, 
as  shown  in  Section  4,  the  revised  estimate  at  the  decision  time  t  is  one  of  the 
following: 


8$'  +C,;1  (E'  Khn  -  0'n)  (8A„-  S  A'n) 


o  A  _ 
^  -n 


(I  +JBn)  85' 


(measurement) 


(2.  16) 


(correction) 
(no  action) 


The  scalar  coefficient  a^  is  computed  from 


=  h‘  K‘ 


E' 


Khn  “ 


2hn 


K  ct>' 
n  Znn 


+  a„ 


(2.  17) 


The  rectangular  matrix  J  has  six  rows  and  three  columns 


(2.  18  ) 


and  is  just  the  reverse  of  the  K  matrix.  The  matrix  Bn  is  also  rectangular 
having  three  rows  and  six  columns  and  is  partitioned  as  shown 


6„=|Cn*  -If  (2.19) 

where  C*  is  one  of  the  fundamental  navigation  matrices  described  in  Section  3.2. 

At  each  decision  point  it  is  also  necessary  to  update  the  correlation 

matrix  E  and  to  compute  the  correlation  vector  0  1  =  <t>  .  0_  _  , . 

n  —  nn  n,  n-l  — n,  n-l 

Thus 


E 


n 


-  %  K  *b,  -  *»)  Kh„ 

+  J  v„  jT 


(measurement) 
(correction) 
(no  action) 


(2.  20  ) 


^n,  n-l 


_  n  -  2 

=  an  -n  -  1  +  £ 

t  =0 


I 

) 


k 


n  -  1 


(2.  21) 
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where 


i  V  (Ek  Khk-0kVbk  xT 


(measurement) 
(correction) 
(no  action) 


(2.  22) 


!al  H  (E£  Kh£  ~±'i0 

]  a  i 

0 


(measurement) 
(correction) 
(no  action) 


(2. 23) 


The  vector  7]  ^  is  the  difference  between  the  commanded  velocity  correction  and 
the  actual  velocity  change  implemented  at  time  t£  . 

The  above  collection  of  formulae  provides  the  means  of  maintaining  an  up 
to  date  estimate  of  the  deviation  vector  6xn  but,  in  themselves,  do  not  provide 
any  clue  as  to  what  decision  should  be  made  at  each  point.  Suggestions  for 
reasonable  decision  rules  are  discussed  in  Section  6.  2  and  in  Appendix  B. 

It  is  important  to  note  that  the  navigation  formulae  are  considerably 
simplified  if  the  measurement  errors  are  considered  uncorrelated.  For  this 
special  case  all  the  0  correlation  vectors  are  identically  zero  so  that  far  less 
computation  is  required. 
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3.  FUNDAMENTAL  NAVIGATION  MATRICES 


Basic  to  the  solution  of  the  navigation  problem  is  a  certain  collection  of 
matrices.  The  objective  here  is  to  introduce  these  matrices,  indicate  their  role 
in  the  navigation  theory,  and  show  how  they  may  be  obtained  as  solutions  of 
differential  equations. 


3.  1  General  Solution  of  the  Linearized  Trajectory  Equations 


Let  rg(t  )  and  v  (t)  denote  the  position  and  velocity  vectors  of  the  space¬ 
craft  in  an  inertial  coordinate  system,  and  let  g  (r  ,t)  denote  the  gravitational  [ 


acceleration  at  position  _rg  and  time  t. 


d-r* 


dt 


-  y. 


Then 


dt 


£  (£*,») 


(3.  1) 


are  the  basic  equations  of  motion  of  the  spaceship  except  for  those  brief  periods 
during  which  propulsion  is  applied. 


Let  the  vectors  _rQ(t)  and  yQ(t)  represent  the  position  and  velocity  at  time 
t  associated  with  the  prescribed  reference  trajectory,  and  define 


sl(0  =_rs(t)  “I0(t)  ,  5y(t)  =vs(t)  -v0(t)  (3>2) 

Then,  the  deviations  6r  and  6v  may  be  approximately  related  by  means  of  the 
linearized  differential  equations: 


d(Sj) 

dt 


8v  ,  -^=G(r0,t)Sr 


(3.  3) 


where  G(r  t)  is  a  matrix  whose  elements  are  the  partial  derivatives  of  the 

components  of  g(r  ,  t)  with  respect  to  the  components  of  r  . 

-o  — o 

A  particularly  useful  fundamental  set  of  solutions  of  Eqs.  (3.  3)  may  be 

developed  in  the  following  way.  Let  t^  and  t^  be,  respectively,  the  time  of 

launch  and  the  time  of  arrival  at  the  target.  Then,  define  the  matrices  R(t ) , 

* 

R  (t)  ,  V(t),  V  (t)as  the  solutions  of  the  matrix  differential  equations 


dR 

dt 


=  V 


(3.  4) 
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which  satisfy  the  initial  conditions 


R (tL)  r  0  ,  R*(»a)=0 

(3.5) 

V(»L)  =  I  ,  V*(»A)=I 

Here  O  and  I  denote,  respectively,  the  zero  and  identity  matrix.  If  we  now  write 

8r(»)  =  R(0  c.  +  R*(t)  c*  (3.6) 

8v(t)  =V(t)  c  +  V*(t)  c*  <3-7> 


where  c  and  c*  are  arbitrary  constant  vectors,  it  follows  that  these  expressions 
satisfy  the  perturbation  differential  equations  (3.  3)  ,  and  contain  precisely  the 
required  number  of  unspecified  constants  to  meet  any  valid  set  of  initial  or  bound¬ 
ary  conditions. 


The  elements  of  the  R  and  V  matrices  represent  deviations  in  position  and 
velocity  from  the  corresponding  reference  quantities  as  the  result  of  certain 
specific  deviations  in  the  launch  velocity  from  its  reference  value.  For  example, 
the  first  columns  of  these  matrices  are  the  vector  deviations  at  time  t  due  to  a 
unit  change  in  the  first  component  of  the  velocity  at  time  t^.  Corresponding 
interpretations  may  be  ascribed  to  the  other  columns  as  well.  A  similar  dis¬ 
cussion  will  provide  a  physical  meaning  for  the  elements  of  R*  and  V*  .  For 
this  purpose,  however,  it  is  convenient  to  imagine  the  roles  of  launch  and  target 
points  as  reversed. 

3.  2  The  Vector  Velocity  Correction 

Associated  with  the  position _rg  and  the  time  t  is  the  vector  velocity  re¬ 
quired  by  the  spacecraft  to  travel  in  free  fall  from_rs(t)  to  the  target  point 
-o  ^  A  ^  *n.  exPress^on  for  this  velocity  vector  is  readily 

obtained  from  Eqs.  (3.  6)  and  (3.  7).  The  condition  that  the  vehicle  pass 
through  the  target  point  is  met  by  the  requirement 

hit  a)  =  0  =  R  (f  A)  fi.  +  R*  (tA)c* 

Since  R*  (t^  )  =  O,  it  follows  that  c  =  0  ,  Eliminating  c_*  between  Eqs.  (3.6) 
and  (3.7)  gives  for  the  required  velocity  deviation*  at  time  t 

8v+(t)  =  V*(t)  R*(t)'’  S^t)  (3,8) 


*The  superscripts-  and  +  are  used  to  distinguish  the  velocity  just  prior  to 
correction  from  the  velocity  immediately  following  the  correction. 
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(3.9) 


* 

Hence,  the  required  velocity  correction  A  v  is  given  by 

Av*(t)  =C*(t)  S_r(t)  -  Sv-(t) 
where  the  C  *  matrix  is  defined  by 

C*(»)  =  V*(t)  R*(t)'’  (3.10) 

The  elements  of  the  C*  matrix  are  deviations  in  vehicle  velocity  from  the 
reference  values,  asrequiredto  place  the  vehicle  on  a  trajectory  to  the  target 
point,  which  arise  from  certain  specific  deviations  in  the  vehicle  position.  The 
interpretation  applied  to  the  columns  is  made  in  the  manner  described  earlier 
in  connection  with  the  R  and  V  matrices. 

If  the  spacecraft  has  been  in  a  free-fall  status  since  launch,  then,  by 
employing  arguments  similar  to  those  used  in  establishing  Eq.  (3.  8)  ,  it  can  be 
shown  that 

Sv-(t)  =C(»)  8r(t)  (3.  11) 

where 

C(t)  =  V(t)  R(t)‘ 1  (3..  12) 

In  this  case  Eq.  (3.9)  takes  the  form 

Av*(t)  =  CC*(t)  -  C(f)3  Sr (t)  (3.  13) 

Since  6r(t)  is  different  from  zero  solely  as  a  result  of  an  injection  velocity 
error  6v(tL  ),  it  follows,  from  the  definition  of  the  R  matrix,  that 

Ay*(f)  =-  A(t)  6v(tL)  (3.14) 

Thus,  the  A  matrix,  defined  by 

A(t)  =V(t)  -  C*(t)  R(t)  (3.15) 

relates  a  deviation  in  launch  velocity  to  the  velocity  impulse  required  at  time  t. 
A  starred  form  of  the  A  matrix 

A*(t)  =  V*(t)  -  C (t)  R*(t)  (3.16) 

will  occur  in  the  subsequent  discussions. 

3.  3  Differential  Equation  Solutions 

The  matrices  C,  C*,  A  ,  A*  may  also  be  generated  directly  as  solutions 
of  differential  equations.  However,  for  C  and  C* ,  a  difficulty  arises  in  pre- 

* 

scribing  appropriate  initial  conditions.  From  the  initial  values  of  the  R  and  R 
matrices,  it  follows  that  C(t^  )and  C*  (t^  )are  both  infinite.  The  singularities 
may  be  avoided  by  working  directly  with  the  differential  equation  for  the  inverse 
matrices  C  *and  C* 
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By  differentiating  the  identity 

v 1 


C(t)  '  V(t)  =R(t) 


,-l 


and  using  Eq.  (3. 4)  ,  the  following  equation  for  C  results 

.-1 


dC  +  C'1  G  C"1 


dt 


Similarly,  we  obtain 


+  C*'1  G  C*'1 


(3.  17) 

(3.  18) 


(3.  19) 


Equations  (3.  18)  and  (3.  19)  may  be  used  to  demonstrate  an  interesting 
property  possessed  by  C  and  C  .  It  is  easy  to  show  that  the  G  matrix  is  sym- 
metrical.  It  follows  at  once  that  the  matrices  C  and  C  will  be  symmetrical  for 
all  values  of  t  in  the  interval  (t^,  t^)  if  they  are  symmetrical  for  any  particular 
time.  But  from  Eq.  (3.  17)  and  a  similar  one  involving  starred  matrices,  we 
have 

C (tL)*  ^  =  0  ,  C*(tA)‘1=0  (3.20) 


so  that  C.  and  C*are,  indeed,  symmetrical  for  t  equal  to  t^  and  t^  respectively. 
Hence  C(t)  and  C*(t)  are  symmetrical  for  all  t  in  the  interval  from  launch  to 
the  target  point. 

In  an  entirely  analogous  manner,  differential  equations  may  be  developed 
for  A  and  A*.  By  differentiating  Eqs.  (3.15)  and  (3.  16)  and  using  Eq.  (3.4), 
one  readily  obtains  the  equations 


~r  + 

dt 


C*A 


=  0 


(3.  21) 


and 


with  the  initial  conditions 


d  A  * 
dt 


+  C 


A*  =0 


A(tL)  =  I  ,  A*(tA)  r  I 


(3.  22) 

(3.  23) 


3.4  The  State  Transition  Matrix 


Let  6r  =  6r(t  )  and  6v  =  Sv(t  )  be  the  deviations  in  position  and 
-  n  -  n  ~n  ~  n 

velocity  at  time  tn,  and  let  Rn,  V  ,  ...  be  the  corresponding  values  of  the 
fundamental  matrices.  The  c  and  c*must  be  obtained  as  solutions  of 


8-rn 


+ 


R 


* 

n 


* 


c 


(3.  24) 


8-n 


=  V„  c 


.  .*  ♦ 
+  Vn  c 


(3.  25) 
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Multiplying  Eq.  (3,24)  by  R~\  we  obtain  for  c 

c  =R;1  (S_r„  -  R*  c*)  (3.26) 

Then,  by  substituting  this  expression  into  Eq.  (3.  25)  and  using  Eqs.  (3.  12)  and 
(3.  16)  ,  there  results 

£*  ~  ”  A  n  (Cn  S-rn  "  (3.27) 

Finally,  from  Eq.  (3.26)  we  have 

A  •  1  /r*  c  j  \  (3.28) 

c  -  -  An  {C„  Srn  -  Syn) 

after  some  simplification.  Thus,  with  c  and  c*  determined,  the  position  and 
velocity  deviations  at  any  other  time  t  are  given  by  Eqs.  (3.  6)  and  (3. 7)  . 

In  terms  of  the  six  dimensional  deviation  vector  defined  by  Eq.  (2.2)  ,  the 
result  may  be  written  in  the  form 

Sx(t)  =  „  . 

-  "  "  (3.29) 

Consider  now  a  specific  value  of  t  =  t  .  Then  substituting  from  Eqs.  (3.  27) 
and  (3.  28)  into  Eq.  (3.  29)  ,  a  relationship  between  6xn+^  and  6xfi  is  displayed 

S?n+1  =  $n+l,n  S-Xn  (3.30) 

where  4>n+j  n,  the  six-dimensional  state  transition  matrix,  is  computed  from 


R(f) 

R*(f) 

c 

V(t) 

V-tt) 

* 

c 

Rn  +  1 

r:+i 

-a;j  o 

C  - 1 

Vn+1 

v:+i 

*-i 

o  -An 

<=n  -  ' 

(3.  31) 
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4.  DERIVATION  OF  THE  OPTIMUM  LINEAR  ESTIMATE 
As  noted  in  the  Introduction,  the  optimum  linear  estimate  of  the  deviation 

A 

vector  may  be  expressed  as  a  recursion  formula.  Therefore,  assume  j 

and  E^  are  known  and  that  a  single  measurement  of  the  type  described  in 

Appendix  A  is  made  at  time  t  .  The  observed  deviation  in  the  measured  quantity 

A  is  6  A  ,  and  the  best  estimate  for  6  A  ,  as  obtained  from  the  extrapolated 
n  n'  A  n'  r 

estimate  of  6xn  is  given  by  Eq.  (2.  15).  Then  a  linear  estimate  for  the 

deviation  vector  6x  at  time  t  is  expressible  as  a  linear  combination  of  the 
~  n  n 

extrapolated  estimate  of  6xn  ^  and  the  difference  between  the  observed  and 
estimated  deviations  in  the  measured  quantity  A  .  Thus 

Sx„  +  Wn  (8A„  -  Va'j  (4.1) 

where  the  vector  wn  is  a  weighting  factor  which  will  be  chosen  so  as  to  minimize 
the  mean-squared  error  in  the  estimate. 

For  this  purpose  use  Eqs.  (2.9),  (2.7)  and  (2.  5)  to  write 
en(wn)  =  Sxn  -  8x„ 

(4.2) 

=  $Xn  +  -n  (SAn  +  an  *  5An>  "  K 
=  d  “  Wn  hj  KT)(8x,;  -  8x„)  +  wn  an 
=  d  “  w„  bl  KT)  e,;  +  w„  an 

where  I  is  the  six-dimensional  identity  matrix.  Then  the  correlation  matrix 

En  defined  by  Eq.  (2.  11)  may  be  expressed  as  a  function  of  the  weighting  vector 

w  as 
—  n 

En(wn)  =d  -wnhnT  KT)  E'  (I  -K  hn  wj)  (4*3) 

+  (I  “  Wn  bj  KT)  0n'n  wj 

+  -n  C  ‘  Kbn  +  -n  an 

The  mean-squared  errors  in  the  estimate  of  position  and  velocity  devia- 
2  2 

tions  and  6^  are  simply  the  respective  traces  of  the  submatrices 
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E(1)  and  E(4). 
n  n 


If  the  six-dimensional  weighting  vector  vv^  is  partitioned  into  two 


three-dimensional  vectors 


(1) 


(2) 


(4.4! 


then  from  Eq.  (4.3)  it  is  easy  to  show  that  E  is  a  function  only  of  w  and 

(4)  (2)  n  n 

E  is  a  function  only  of  w  n  .  Therefore,  for  the  purposes  of  the  following 

discussion,  it  is  legitimate  formally  to  treat  the  mean- squared  error  in  the 
9 

estimate  e  (w  )  as  the  trace  of  the  six-dimensional  correlation  matrix  E  (w  ) 
n  — n  n  — n 

The  subvectors  of  the  optimum  weighting  vector  w^  will  then  each  be  optimum 
for  the  respective  estimates  of  position  and  velocity  deviations. 

In  order  to  determine  the  optimum  weighting  vector,  one  may  apply  the 
usual  technique  of  the  variational  calculus.  Let  take  on  a  variation  6wn  and 


obtain  from  Eq.  (4.3) 

«<!< *„>  =  2  »  [  -  s»„b„TKTE^  (I  -Kb„«J) 


(4.  5) 


'  8!?n  hi  K^nn  *1  +  S*n  £n l  (•  *  K  K  S J>  +  «!?„  *1  an  ] 


If  6  e  (w  )  is  to  vanish  for  all  variations  6  w  , 
-  n  — n 


n  r.n  ‘ 

then  it  must  follow  that 


°n  -n  =  En  K  bn  "  £nn 

where  the  positive  scalar  quantity  an  is  defined  by  Eq.  (2.17) 


(4.  6) 


It  can  be  readily  shown  that  the  \yn  determined  from  Eq.  (4.6)  actually 
does  minimize  en^n^  •  Suppose  that  the  optimum  is  replaced  by  another 
weighting  factor  w^  -  y^ 


Then  from  Eqs,  (4.3)  and  (2.17) 

J( w„  -  Yt)  =fr  [En  -  2  (wn  -  yn)  (hj  KT  +  an  (w„  -_yn)  (wj  -  yj)] 

and  using  Eq.  (4.  6) 


(4.  7) 


so  that 


®n  *  In)  =  tr  [  En  "  an  "  vJ  1*1  +  vj)] 
en  "  Vn)  =  en  (*„>  +  °n  tr  tyn  vj) 


(4.  8) 


(4.9) 

Thus,  the  mean-squared  error  is  not  decreased  by  perturbing  wn  if  Eq.  (4.  6) 
holds. 


Having  obtained  the  optimum  weighting  vector,  the  expression  for  the 
correlation  matrix  of  the  estimate  errors  En  given  by  Eq.  (4.  3)  may  be  written 
in  a  more  convenient  form.  Thus,  from  the  definition  of  a^  in  Eq.  (2,17), 
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there  results 

E„  =  En  n  -  K  hn  wj)  -  wn  hi  KT  E '  +  wn  0'J  +  <p'  wl  +  on  wn  wl 

n  n  - n  - n 1  -n  n  n  -n  inn  _ nn  — n  n  -n  —  n 

Substituting  from  Eq.  (4.6),  the  final  expression  may  be  written  as 

e„  (e;  Khn  -  *;„)  (e;  Kh„  -*yT 


(4.  10) 
(4.  11) 


For  the  case  in  which  the  measurement  errors  are  uncorrelated,  i.  e. , 

0'  =0  for  all  n  and  m,  the  estimation  procedure  is  complete.  Equations 

”  rim 

(4.  1)  and  (4.11)  in  the  form 

*«„  =  «*„'  *»;’  K  <SA„  -  «A')  <4-12) 

E„  =E;  <E'  K b„) (E •  Kb„)T 

then  serve  as  recursive  relations  to  be  used  in  obtaining  improved  estimates 
of  position  and  velocity  deviations  at  each  of  the  measurement  times  t^,  t^,  .  .  . 


If  the  measurement  errors  are  correlated,  the  procedure  must  be  expanded 
to  include  a  method  for  computing  0*  .  For  this  purpose,  a  recursion  formula 

may  be  developed  from  Eq.  (4.2)  .  We  have  at  time  t  ^ 


-em  - 1  =  0  -  - 1  hi .  ]  kt)  e; . ,  + 


(4. 13) 


so  that  j  T  _ 

^nm  m,m  *  1  ^  -m-1  -m  -1  ^  )  ^n,m  -  1  +  ^  m,m  -  1  -m-1  an  am -  1 


(4. 14) 


m  =  2,  3,  ...»  n 

Equation  (4.  14)  permits  the  calculation  of  0  1  by  successive  substitu¬ 
tion  beginning  with  m  =  2  and  noting  that 


^1,1  =<1>1,0  an  ®c  (4.15) 

where  e  =  -  is  the  negative  of  the  deviation  vector  at  time  t2  arising  from 
improper  injection  into  orbit. 


The  closed  form  calculation  given  in  Eq.  (2.  21)  is  readily  derived  from 
Eq.  (4.14)  . 
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5.  STATISTICAL  ANALYSIS  OF  THE  GUIDANCE  PROCEDURE 


From  exact  knowledge  of  the  six-dimensional  deviation  vector  6  xn  at  time 
t  ,  a  velocity  correction  may  be  calculated  which,  if  implemented,  will  insure 
the  vehicle's  arrival  at  a  fixed  point  in  space  at  the  required  time.  However, 
only  the  estimate  6x^  is  available.  From  this, an  estimate  of  the  velocity  corr¬ 
ection  vector  A  vn  may  be  determined  from 

Avn=BnSin  (5.1) 

where  is  defined  by  Eq.  (2.  19)  .  (Refer  to  the  discussion  leading  to  Eq.(3.9). ) 

The  need  for  a  velocity  correction  arises  solely  from  improper  injection 
into  orbit.  If  the  first  such  correction  is  executed  perfectly,  then,  of  course, 
no  further  corrections  are  required.  However,  because  of  imperfect  knowledge 
of  position  and  velocity  obtained  from  navigational  measurements,  the  command¬ 
ed  velocity  change  will  be  in  error.  Furthermore,  the  actual  velocity  change 
experienced  will  differ  from  that  commanded  because  of  imperfect  instrument¬ 
ation.  Therefore,  subsequent  corrections  will  be  required  to  remove  the 
effects  produced  by  earlier  inaccuracies. 

5. 1  Correlation  Matrix  of  the  Velocity  Correction  Vector 


For  notational  purposes  in  this  section,  it  will  be  convenient  to  distin¬ 
guish  the  times  of  velocity  correction  from  the  other  decision  times.  Therefore, 

t  will  be  used  to  denote  the  time  of  the  n-th  corrective  maneuver.  An  arbi- 
c,  n 

trary  number  of  decision  points  may,  of  course,  fall  between  t  and  t 

c,  n  c,  n+ 1 

Suppose  that  corrective  action  is  to  be  taken  at  the  time  t  .  Let  the 

A  ^ 

commanded  velocity  change  be  A  yc  n>  while  the  actual  velocity  change  exper¬ 
ienced  is  Av  .  Then,  with  n  denoting  the  uncertainty,  we  have 
Gj  ii  v  c»  n 


aO 

-c,n 


=  Avc  _ 

-c,n 


'e,n 


(5.2) 


The  actual  velocity  change  may  be  expressed  as 

A  v  _  _  =  B_  n  (Sx'_  +  e '  -  -n 

_ c,n  c, n  '  _  c,n  —  c^n7  j  cfn 


(5.  3) 
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where  e'c  is  the  error  vector  as  extrapolated  from  the  last  observation  point. 
The  minus  superscript  is  used  to  emphasize  that  the  deviation  vector  corres¬ 
ponds  to  the  time  immediately  prior  to  application  of  the  corrective  velocity 
impulse.  Similarly,  a  plus  superscript  will  distinguish  a  deviation  vector 
immediately  subsequent  to  a  velocity  correction. 


Just  as  the  extrapolated  error  vector  is  altered  at  an  observation  point, 
so  also  will  it  change  at  a  correction  point.  Thus 

(5.4) 


6  —  0  + 
-c,n  =  c,n 


-  c,n 


so  that  Eq.  (5.3)  becomes 

A*c,n  =  Bc,n  (8*c,n  +  Sc,r»)  (5.5) 

We  further  note  that  the  correlation  matrix  of  the  deviation  errors  must  be 
updated  at  a  correction  point  following  the  application  of  the  velocity  correction. 
It  follows  from  Eq.  (5.4)  that 

(5.6) 

II  II 

E.,n  =E;.n  * 


0 

0 

0 

T 

3c, n  ?c,n 

A  more  convenient  form  of  Eq.  (5.5)  is  necessary  for  subsequent  statis¬ 
tical  analysis.  Indeed,  it  is  possible  to  express  the  velocity  correction  actually 

applied  at  time  t  directly  in  terms  of  the  error  vectors  e  and  e 

c,  n  C|  n  Cj  n - 1 

For  this  purpose,  we  use  the  definition  of  the  transition  matrix  operator  to 
write  Eq.  (5.  5)  as 

A-e,n  ~  Bc,n  ^c,n;c,n-l  8-c,n-l  +  —  c,n^ 


Bc,n  ^c,n;n  -  I  8-c,n  -  l  + 


A¥c,n-1 


>+  B. _  e.  „ 

Cj  n  C|  n 


(5.  7  ) 
(5.8) 


Further  simplification  is  possible  because  of  the  identity 

B  $  i  =  /\  /\  i  5  i 

uc,n  c,rt,c,n-l  c,n  c,n  -  1  c,n-I 
which  is  readily  established  from  Eq.  (3.  31)  . 

Therefore,  .  i 

A*c,n  =  Ac,n  Ac,n-1  (Bc,n-1  8 *'c,n  -  1  '  A*c,n  -  l)  +  Bc,n  ®c,n  (5-9) 
However,  since  Eq.  (5.  5  )  obtains  at  the  previous  correction  time  t£  n_j.  then 

Bc,n  -  I  8*c,n-1  "  8-c,n  •  1  =  "  Bc,n  -  1  ?c,n  -  1  (5.  10) 
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Hence,  finally,  we  have 


-  1 


Av„  _  =B__e,_  -  A  A’’  i  Br  i  e  i 

c,n  c^n  I  c,n  -  I  n  * 


J-c,n  “c,n-c,n  '  '  c,n  '  'c,n  -  I  “c,n  -  I  —  c,n-l  (5.11) 

as  the  desired  expression  or,  alternately,  using  Eq.  (5.  2)  ,  we  may  write 


^~c,n  Bc,n  -c,n  ^c.n^c.n-l^c.n-l-c.n-l  (5.12) 

The  correlation  matrix  of  the  estimated  velocity  correction  vector  is 
found  by  computing  the  mathematical  expectation  of  the  product  of  A  vc 
and  its  transpose.  Thus 


*?«..  <»  -  ®c,n  ®c,n  -J,n  •  1  ®e,n-  ,  ) ,  aJ,„ 


A,  _  A 


-  1 


<T  dT 


c,n  /'c,n-l  Dc,n-1  -c,n-l  -c,n  Dc,n 


+  Ac,n  Ac'n-  1  BC,n-l  Ec,n-1  BI,n  - 1  Ac,n  -  1  Ac,n 


(5. 13) 


expresses  this  correlation  matrix  in  terms  of  the  correlation  matrices  of  the 

errors  in  the  deviation  vectors  at  the  two  consecutive  correction  points  and  a 

T 

cross-correlation  matrix  e1  e 

-  c,  n  -  c,  n-  1 

The  cross-correlation  matrix  appearing  in  Eq.  (5.  13)  is  most  easily- 
computed  by  a  recursive  operation.  In  order  to  facilitate  a  description  of  the 
procedure,  let  us  identify  the  time  t^  with  fc  n  ^  and  assume  p  observation 
times  before  the  next  velocity  correction  at  time  t^  Then,  from  the  defin¬ 
itions  of  the  quantities  involved. 


^m+l-m  +  Em 


(5.  14) 


At  the  k-th  observation  point,  the  error  vector  associated  with  the  deviation 
vector  is 

fim+lc  =d  -  +lc  i>m  +  k  kT)  +  k  +  «m+k  am  +  k  (5.15) 

according  to  Eq.  (4.2)  .  Therefore,  in  general,  we  have 

>  T  f  /*  i  T  i,  T  \  ,  T  ( 5 .  16) 

—  m  +  k  +  1  -m  _<^m+k+l,m+k  ^m+k-m+k^  y^m+k-m  +  -m+k-m+k,m 

Now,  starting  with  Eq.  (5.  14)  ,  a  succession  of  cross-correlation  matrices  is 

p.  The  final  calculation  produces 
T 


computed  from  Eq.  (5.  16)  for  k=  1,2, 

77  ~T~ 


-  m+p+I  -  m 
the  evaluation  of  Eq.  (5.  13) 


which  is  identical  with  the  matrix  e 


—  c,  n  —  c,  n-  1 


required  for 
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The  mean- squared  estimate  of  the  velocity  correction  is  determined  as 
A  AT 

the  trace  of  the  matrix  Av  A  v  .  As  a  basis  for  a  decision  theory,  it  is 

-  c,  n  c,  n 

important  to  know  something  of  the  precision  of  the  estimate.  Clearly,  a 
velocity  correction  having  a  large  uncertainty  should  not  be  commanded  if  it  is 
possible  to  improve  substantially  the  estimate  by  future  observations.  The  un- 

A 

certainty  cjc  ^  in  the  estimate  Av^  ^  is  simply 


-c,n 


A?c,n  "  Bc,n  k-c,n  =  Bc  n  e'f„ 


(5. 17) 

Hence,  the  mean-squared  uncertainty  is  determined  as  the  trace  of  the  matrix 


5.  2 


d  dT  E'  BT 

-c,n-c,n  c,n  ,_c,n  c,n 


Uncertainty  in  the  Applied  Velocity  Correction 


(5.  18) 


In  order  to  complete  the  statistical  analysis  of  the  velocity  correction,  it 
is  necessary  to  examine  more  carefully  the  vector  uncertainty  rj  in  the  velocity 

—  A 

correction.  The  inaccuracy  in  establishing  a  commanded  velocity  correction  Av 
is  due  to  errors  in  both  magnitude  and  orientation.  In  the  following  analysis  the 
two  sources  of  error  will  be  assumed  independently  random  with  zero  means. 


Consider  a  coordinate  system  in  which  the  estimated  velocity  correction 
vector  is  along  one  of  the  coordinate  axes.  Then  if  M  is  the  transformation 
matrix  which  relates  the  selected  axis  system  and  the  original  reference  system, 
we  may  write 

Ay  ~  Ay  M 

Now,  define  a  random  variable  k  such  that 

Ay  =  (1  +  k 

and  let  y  be  the  random  angle  between  A  v  and  Ay  It  will  be  assumed  that  both 
k  and  y  are  small  quantities  so  that  powers  and  products  are  negligible  compared 
with  unity.  The  actual  vector  velocity  correction  is  then 

Av  =  (1  +  k)  Av  M 


y  cos  P 
y  sin  P 

1 


)  Av 


(5.  19) 


(5.  20  ) 


where  P  is  a  polar  angle  defining  the  rotation  of  Ay  with  respect  to  A  v  .  Hence, 
the  uncertainty  vector  rj  is  expressible  as 


Av 


Av 


-  Ay  M  |  (1  +  k)  y 


cos  p 

0 

sin  P 

+  K 

0 

0 

1 

(5.22) 
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Assume  that  * ,  j  ,  {3  are  statistically  independent  random  variables  with 
zero  means.  Further  assume  that  (3  is  uniformly  distributed  over  the  interval 
0  to  2  .  Then  one  obtains  for  the  correlation  matrix  of  the  velocity  correction 

uncertainty 


vv 


k  Av  Ay  +  Av^  M 


1 

0 

0 

0 

1 

0 

0 

0 

0 

(5.  23) 


_  2  aa 

—  K  Av 


.  aT  y  /,aT  .a  .  a  .  aT  \ 
Ay  +  — y  (^Ay  Ay  I  -  Ay  Ay  J 


whei  a  I  is  the  three-dimensional  identity  matrix  and  *  and  }  are  the  mean- 
squared  valued  of  k  and  y . 


5.  3  Miss  Distance  at  the  Target 


Turning  now  to  the  problem  of  guidance  accuracy,  the  determination  of 
the  position  deviation  vector  at  the  nominal  time  of  arrival  at  the  target  is  made 
by  extrapolating  the  deviation  vector  from  the  point  of  the  final  velocity  corr¬ 
ection.  Thus,  if  t^  is  the  time  of  the  last  correction  and  6  x.^  is  the  deviation 
vector  at  the  time  of  arrival  t^,  then 


S*A  '  ^A,N  S*N 


(5.  24) 


But  from  Eq.  (3.  31)  and  the  terminal  conditions  for  the  navigation  matrices, 
we  have 


^A,N 


“  R  A  A 


-  1 
N 


_VA  An' 


-  A 


*  - 1 


r* 


(5. 25) 


Hence,  the  position  deviation  vector  at  the  target  6_r^  may  be  written  as 


slA  ~  "RA  AN  BN  5*N 


(5.  26) 


with  a  similar  expression  obtainable  for  the  velocity  deviation  at  time  t^. 

The  target  position  error  may  be  written  ultimately  in  terms  of  the 
error  vector  _e^T  according  to  the  following  self-evident  steps 
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SlA  "  -  RA  AN  BN  <ix~N  T  J  -YN) 


■  -  Ra  AN  (Bn  *xn  -  AvN) 


RA  aN  <Bn  5N> 


Ra  an  Bn  eN 


(5. 27) 


The  mean  square  position  error  at  the  target  is  then  computed  as  the  trace  of 

the  matrix  6r»  6  r  T. 

A 
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6.  APPLICATION  TO  TRANS-LUNAR  NAVIGATION 


6.  1  Decision  Rules 

As  a  necessary  step  in  the  application  of  the  navigation  and  guidance 
scheme  formulated  in  this  paper,  certain  rules  must  be  adopted  concerning  the 
course  of  action  to  be  taken  at  each  of  the  "decision  points"  described  in 
Section  2.3.  The  number  and  frequency  of  observations  must  be  controlled  in 
some  manner  --  ideally  by  a  decision  rule  which  is  realistically  compatible  with 
both  the  mission  objectives  and  the  capabilities  of  the  measuring  device.  If  an 
observation  is  to  be  made,  a  decision  is  required  regarding  the  type  of  measure¬ 
ment  and  the  celestial  objects  to  be  used.  Periodic  velocity  corrections  must 
be  applied  and  the  number  of  impulses  and  times  of  occurrence  must  be  decided. 

Once  the  decision  rules  have  been  specified,  it  is  necessary  to  test  their 
effectiveness  according  to  some  measure  of  performance.  A  typical  objective 
is  to  minimize  the  miss  distance  at  the  target.  However,  a  reduction  in  miss 
distance  usually  implies  an  increase  in  either  the  required  number  of  measure¬ 
ments  or  a  greater  expenditure  of  corrective  propulsion  or  both.  In  the  face  of 
these  conflicting  objectives,  compromises  are  clearly  necessary  and  statistical 
simulation  provides  a  means  of  arriving  at  an  acceptable  balance. 

In  the  interest  of  minimizing  the  number  of  simulator  runs,  Monte  Carlo 
techniques  should  be  avoided  if  possible.  Fortunately,  it  is  unnecessary  to 
generate  the  true  spacecraft  trajectory,  as  would  be  required  for  Monte  Carlo 
simulation,  in  order  to  analyze  the  effects  of  a  particular  set  of  decision  rules. 
The  reader  may  readily  verify  that  Eq.  (2.  16) ,  which  defines  the  estimate 
6  xn  and  depends  on  actual  measurement  data,  is  never  involved  in  any  of  the 
statistical  calculations. 

A  specific  example  of  a  set  of  decision  rules  to  be  applied  at  each  decision 
point  is  as  follows:  _ 

1.  The  estimated  mean-squared  velocity  correction  Av^  and  the  mean- 

2  n 
squared  uncertainty  d^  associated  with  the  estimate  are  computed  from 

Eqs.  (5.  13  )  and  (5.  18).  If  the  ratio 


(6.  1) 


is  less  than  a  specified  amount  R  ,  ...  a  velocity  correction  is  made  at 

r  v\min) 

time  t  .  In  any  case,  the  correction  is  made,  regardless  of  the  value  of  R^, 

if  the  root-  mean-squared  velocit\’  correction  exceeds  a  fixed  amount  Av 

M  “  max 

2.  If  neither  of  the  criteria  is  met  which  would  call  for  initiation  of  a 

velocity  correction,  the  desirability  of  making  an  observation  is  examined.  For 

this  purpose,  an  abbreviated  star  catalog  is  postulated  together  w  ith  selected 

planets.  Each  star  and  planet  combination  is  analyzed  to  determine  the  effect 

of  each  such  measurement  on  the  potential  reduction  in  the  miss  distance  at  the 

target.  More  specifically,  at  time  t  the  mean- square  miss  distance,  which 

would  result  if  no  further  corrective  action  were  taken,  is  calculated  from 


Eq.  (5.27)  .  Now,  if  any  measurement  were  made  and  were  followed  imme¬ 
diately  by  a  velocity  correction,  a  calculable  mean-square  reduction  in  the 
target  error  would  result.  The  particular  star-planet  combination  producing 
the  greatest  mean- square  reduction  in  target  error  is  then  defined  as  the 
best  potential  measurement. 

2  +  2- 

Now  let  6r^  and  6r^  be  the  respective  mean-square  miss  distances 
which  would  result  with  and  without  the  best  possible  observation.  Then,  if  the 

(6.2) 


is  greater  than  a  specified  value  Rp(max)>  t^ie  best  potential  measurement  is 

made  at  time  t  .  In  other  words,  for  a  measurement  to  be  made,  a  significant 

reduction  in  the  potential  miss  distance  must  result.  If,  on  the  other  hand, 

the  above  criterion  is  not  met,  no  action  is  taken  at  the  decision  point  t  . 

n 

6.  2  A  Numerical  Example 


The  particular  set  of  decision  rules  formulated  in  the  previous  section 
were  applied  to  the  problem  of  navigating  along  a  trans-lunar  trajectory.  For 
simplicity,  possible  cross- correlation  between  measurement  errors  was 
ignored.  Furthermore,  only  the  Earth  and  the  Moon,  together  with  the  ten 
brightest  stars,  were  considered  for  potential  measurements. 
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For  the  specific  trajectory  chosen  for  illustration,  the  date  and  time  of 
orbital  injection  was  Julian  Day  2440042.  5  at  the  twelfth  hour.  The  target  point 
was  selected  approximately  11,  000  miles  beyond  the  moon  with  the  closest 
point  of  approach  some  3,  000  miles  from  the  lunar  surface.  The  nominal  time 
of  flight  from  injection  was  84  hours. 

The  correlation  matrix  of  injection  errors  Eq  was  assumed  to  be 


which  corresponds  to  assuming  zero  mean-squared  position  error  and  an 
isotropic  velocity  error  whose  root-mean- square  value  is  20  miles  per  hour. 

At  each  decision  point,  twenty  potential  measurements  were  examined 
and  evaluated  according  to  the  decision  criterion.  If  the  best  potential  measure¬ 
ment  was  actually  made,  then  the  second  best  measurement  was  also  evaluated 
using  the  same  criterion.  In  this  manner,  the  possibility  of  two  observations 
at  each  decision  point  was  allowed. 

For  simplicity,  only  star  elevations  above  an  illuminated  horizon  of 
either  the  Earth  or  Moon  were  considered.  Certain  practical  constraints  were 
imposed  so  that  physically  unrealizable  measurements  were  screened  out.  For 
example,  no  measurement  could  be  made  if  the  line  of  sight  to  either  star  or 
planet  edge  were  closer  than  fifteen  degrees  from  the  direction  to  the  Sun. 
Furthermore,  the  relative  orientation  of  the  Earth  and  Moon  W'ere  taken  fully 
into  account;  e.  g,  ,  if  the  illuminated  face  of  the  Moon  formed  the  background 
of  the  edge  of  the  Earth  from  which  a  star  elevation  was  to  be  reckoned,  that 
particular  measurement  would  not  be  made. 


The  optical  measuring  device  used  for  the  observations  was  assumed  to 
be  unbiased  with  a  random  error  whose  variance  was 

cr2  =  (0.00005)2  +(— ■ -f  radians 
v  rSPy 

where  rgp  is  the  distance  in  miles  from  the  spacecraft  to  the  planet  edge.  In 
this  manner  it  was  possible  to  account  for  the  larger  uncertainty  in  defining  the 
horizon  which  would  exist  when  the  spacecraft  is  close  to  the  planet.  At  large 
distances  the  rms  error  is  approximately  0.  05  milliradians. 

Finally,  in  this  example  calculation,  the  error  in  applying  a  velocity 
correction  was  assumed  to  be  isotropic  and  proportional  to  the  commanded 
correction.  Specifically,  the  relation 

=0.0001 

was  adopted  so  that  the  rms  error  would  be  one  percent  of  the  rms  correction. 

Preliminary  results  of  an  analysis  of  this  sample  trajectory  are  summar¬ 
ized  in  the  accompanying  tables.  A  number  of  simulated  guidance  flights  were 
made  for  which  the  strategy  parameters  Rv  and  Rp  had  various  assigned  values. 
Certain  pertinent  navigation  data  are  recorded  in  the  tables  as  functions  of  these 
strategy  parameters. 

In  Table  1  the  navigation  data  is  given  as  a  function  of  the  miss  distance 
reduction  ratio  Rp  for  a  value  of  the  velocity  correction  uncertainty  ratio 

=  0,3.  As  one  requires  each  measurement  to  have  a  proportionately  greater 
significance  in  the  reduction  of  the  potential  target  error,  the  total  number  of 
measurements  is  considerably  lessened.  The  extremes  are  126  and  22  mea¬ 
surements.  The  number  of  corrective  velocity  impulses  remains  about  the 
same  but  the  total  of  the  velocity  corrections  applied  increases  by  fifty  percent. 
The  miss  distance  at  the  target  varies  between  one  and  three  miles  while  the 
uncertainty  in  the  knowledge  of  the  vehicle  velocity  at  the  target  falls  between 
one-quarter  and  one-half  mile  per  hour. 

The  navigation  data  recorded  in  Table  2  is  a  function  of  the  velocity  corr¬ 
ection  uncertainty  ratio  with  Rp  fixed  at  0.4.  The  total  of  77  measurements 
did  not  vary  with  R^  but  the  number  of  velocity  corrections  decreased  from  9 
to  4  as  R  was  decreased  from  0 . 4  to  0  .  1.  The  final  miss  distance  and  uncer¬ 
tainty  in  velocity  were  not  greatly  effected  and  the  total  of  the  velocity  correct¬ 
ions  varied  between  approximately  65  and  80  miles  per  hour.  Also  recorded 
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in  the  table  is  the  time  from  injection  to  the  application  of  the  first  velocity 
correction.  As  the  parameter  R^  was  decreased,  this  time  increased  from  two 
to  five  hours. 

Table  I.  Navigation  data  as  a  function  of  miss  distance  reduction  ratio. 

Velocity  Correction  Uncertainty  Ratio  =  0.3 


Miss 

Distance 

Reduction 

Ratio 

Number 

of 

Measurements 

Number 

of 

Velocity 

Corrections 

Total 

Velocity 

Correction 

(mph) 

Final 

Miss 

Distance 

(miles) 

Final 

Velocity 

Uncertainty 

(mph) 

0.3 

126 

7 

55.5 

1.15 

0.23 

0.4 

77 

7 

71.1 

1.41 

0.28 

0.5 

50 

7 

73.1 

1.90 

0.39 

0.6 

33 

6 

79.7 

2.25 

0.4^ 

0.7 

22 

6 

84.7 

3-19 

0.57 

Table  2.  Navigation  data  as  a  function  of  velocity  correction  uncertainty  ratio. 
Miss  Distance  Reduction  Ratio  =  0.4  No.  of  Measurements  =  77 


Velocity 

Correction 

Uncertainty 

Ratio 

— 

Time 
of  Fi rst 
Velocity 
Correction 

Number 

of 

Velocity 

Corrections 

Total 

Velocity 

Correction 

(mph) 

Final 

Miss 

Distance 

(miles) 

Final 

Velocity 

Uncertainty 

(mph) 

0.40 

2.0 

9 

70.2 

1.41 

0.29 

0.30 

2.0 

7 

71.1 

1.41 

0.28 

0.25 

2.4 

6 

70.6 

1.45 

0.29 

0.20 

2.6 

6 

65.9 

1.37 

0.28 

0.15 

3.4 

5 

74.7 

1.47 

0.30 

0.10 

5.0 

4 

79.1 

1.52 

0.31 

In  order  to  evaluate  the  effect  on  the  navigation  data  of  a  variation  in  the 

time  of  year,  a  number  of  pseudo-trajectories  were  generated  by  the  simple 

device  of  rotating  the  direction  of  the  Sun  as  viewed  from  the  Earth.  The 

trajectory  was  considered  to  be  unchanged  by  this  process  --  the  assumption 

being  quite  adequate  for  the  purpose  of  this  preliminary  analysis.  One  set  of 

values  for  R  and  R  was  selected  and  the  Sun  direction  was  altered  in  sixty 
v  p  J 

degree  steps.  In  this  manner  differ ent  illuminated  portions  of  the  Earth  and 
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Mars  were  visible  to  the  spacecraft  resulting,  thereby,  in  different  measure¬ 
ments.  As  seen  in  Table  3,  the  total  number  of  measurements  varied  by  twenty 
percent  while  the  total  of  the  velocity  corrections  had  a  variation  between  67  and 
87  miles  per  hour.  The  effects  on  miss  distance  and  velocity  uncertainty  were 
not  substantial. 


Table  3.  Navigation  data  for  pseudo-trajectories  as  a  function  of  sun  direction  rotation. 


1  °.5 

Time 

=  Start  to  10  hours 

Miss  Distance  Reduction  Ratio  =  /  q.4 

Time 

=  10  hours  to  50  hours 

l  0.5 

Time 

=  50  hours  to  end 

Velocity  Correction  Uncertainty  Ratio  =0.3 


Sun 

Direction 

Rotation 

(degrees) 

Number 

of 

Measurements 

Number 

of 

Velocity 

Corrections 

Total 

Velocity 

Correction 

(mph) 

Final 

Miss 

Distance 

(miles) 

Final 

Velocity 

Uncertainty 

(mph) 

0 

57 

7 

73.5 

1.70 

■■ 

60 

59 

7 

68.4 

1.82 

120 

53 

6 

67.2 

1.75 

— 

180 

53 

7 

75.8 

2.06 

mm 

240 

55 

7 

75.5 

2.11 

3  00 

49 

6 

86.9 

2.21 

MEM 

Finally,  in  Table  4,  a  complete  history  of  one  of  the  guided  translunar 
flights  is  recorded.  With  =  0.  7  and  Rv  =  0.  3,  the  time  and  nature  of 
each  observation  and  velocity  correction  is  listed.  It  is  curious  to  note  the 
frequency  with  which  a  particular  star  and  planet  combination  is  repeated 
in  the  record  of  observations. 

In  conclusion,  it  should  be  emphasized  that  only  fragmentary  results 
have  been  obtained  by  utilizing  this  navigational  scheme.  Therefore,  too  great 
an  importance  should  not  be  attached  to  the  numerical  data  presented  in  this 
section.  At  this  point,  only  the  ideas  are  significant,  and  one  must  await  a 
more  thorough  investigation  before  definitive  conclusions  can  be  drawn. 
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Table  4.  Example  of  translunar  navigation. 


Miss  Distance  Reduction  Ratio  -  0.7 
Velocity  Correction  Uncertainty  Ratio  =  0.3 


TIME 

(hours) 

OBSERVATION 

VELOCITY  CORRECTION 
(tnph) 

■SB 

Procyon 

■Elf 

Earth 

Procyon 

III 

Earth 

Procyon 

37.05 

mm 

Earth 

Procyon 

3.80 

Earth 

Procyon 

6.00 

Earth 

Procyon 

6.00 

6.80 

Moon 

Capella 

17.14 

.O.'JO 

F.ar.'li 

Procyon 

11.00 

Moon 

Capella 

15.00 

Moon 

Procyon 

18.00 

20.00 

Moon 

Procyon 

6.57 

22.00 

Moon 

Procyon 

28.00 

Moon 

Procyon 

34.00 

Moon 

Procyon 

42.00 

Moon 

Capella 

42.00 

45.00 

Earth 

Arcturus 

3.44 

51.00 

Moon 

Capella 

65.00 

Moon 

Capella 

76.00 

Moon 

Rigel 

80.20 

80.40 

Moon 

Arcturus 

8.31 

81.20 

Moon 

Capella 

82.50 

Moon 

Capella 

12.19 
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APPENDIX  A 


NAVIGATIONAL  MEASUREMENTS 

The  mathematical  processes  are  considered  here  in  some  detail  for 
determining  spacecraft  position  by  means  of  both  celestial  observation  and 
ground  based  radar  measurements.  It  is  assumed  throughout  the  analysis  that 
approximations  to  spacecraft  position  and  velocity  are  already  known  so  that 
perturbation  techniques  may  be  employed. 

Secondary  effects  arising  from  the  finite  speed  of  light,  the  finite  distance 
or  stars,  etc.  are  ignored  in  this  analysis.  Such  effects  may  be  lumped  to¬ 
gether  for  a  particular  reference  point  on  the  trajectory  as  a  modification  to 
the  stored  data  which  represent  reference  values  for  the  quantities  to  be  meas¬ 
ured  at  that  point. 

For  simplicity  in  the  present  analysis,  it  will  be  assumed  that  the 
spacecraft  clock  is  perfect  so  that  all  measurements  are  made  at  known  instants 
of  time.  Methods  of  including  clock  errors  in  the  computation  are  discussed 
thoroughly  in  reference  2. 

As  indicated  in  Section  2.  1  each  measurement  establishes  a  component  of 
spacecraft  position  along  some  direction  in  space.  If  A  is  the  quantity  to  be 
measured  and  6 A  is  the  difference  between  the  true  and  the  reference  values, 
then  it  will  be  shown  that  the  relation  between  6  A  and  the  deviation  in  space¬ 
craft  position  Sr  is 

SA  =hT  Sr  (A  .1) 

regardless  of  the  type  of  measurement.  Thus,  the  ji  vector  alone  will  charact¬ 
erize  the  kind  of  measurement. 

Sun-Planet  Measurement 

The  first  type  of  measurement  to  be  considered  is  that  of  the  angle  from 
the  Sun  to  a  planet.  By  passing  to  the  limit  of  infinite  distance  from  one  or  the 
other  of  these  bodies,  corresponding  relations  for  the  Sun-star  or  planet-star 
type  of  measurement  may  be  obtained. 
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Let  S  and  P  be,  respectively,  the  reference  positions  of  the  space¬ 
craft  and  a  planet  at  the  time  of  the  measurement.  Let  rg  be  the  vector  from 
the  Sun  to  Sq  and  z  the  vector  from  Sq  to  P  .  With  A  denoting  the  angle  from 
the  Sun  line  to  the  planet  line,  we  have 


cos  A  -  -  (r  •  z)/rz  (A.  2) 

where  r  and  z  denote  magnitudes  of  the  respective  vectors  and_z.  Treating 
all  changes  as  first-order  differentials,  it  can  be  shown  that 


8A 


< 


m  —  (n  *  m)  n  n  -  (n  •  m)  m  ' 


r  sin  A 


z  sin  A 


Sr 


(A.  3) 


For  details  the  reader  is  referred  to  reference  2.  Here  m  and  ji  are,  respect¬ 
ively,  the  unit  vectors  from  S  toward  the  Sun  and  toward  P  .  The  two  indivi- 
dual  vector  coefficients  of  6r  in  Eq.  (A.  3)  are  vectors  in  the  plane  of  the  meas¬ 
urement  and  normal,  respectively,  to  the  lines-of-sight  to  the  Sun  and  to  the 
planet . 


Planet  Diameter  Measurement 


If  D  is  the  actual  diameter  of  a  planet,  the  apparent  angular  diameter  A 
is  found  from 


sin  (A/2)  =  D/2z 

Again  taking  differentials  as  before,  one  can  show  that 

Dm  •  Sr 


S  A 


z^  cos  (A/2) 


(A.  4) 

(A.  5 ) 


Star  Occultations 


The  next  type  of  measurement  to  be  considered  is  that  of  noting  the  time 

at  which  a  star  is  occulted  by  a  planet.  Let  z  be  the  vector  from  S  to  P  ,  r  the 

°  o  “ 

vector  from  the  Sun  to  S  and  n  a  unit  vector  in  the  direction  of  the  star  to  be 

o  — 

occulted.  With  y  denoting  the  angle  from  the  star  line  to  the  planet  line  as 
shown  in  Fig.  A-l,  we  have,  at  the  nominal  instant  of  occultation, 

n  •  z  =  z  cos  y  (A.  6) 

Treating  changes  as  first  order  differentials  we  obtain 

■  8z  =  cos  y  Sz  -  z  sin  y  Sy  (A. 7) 


-  cos  y  m  •  8z-zsiny8y 

where  m  is  a  unit  vector  from  S  toward  P  . 

—  o  o 
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Fig.  A-1.  Measurement  of  time  of  a  star  occupation. 


The  angle  deviation  6y  is  computed  from  a  first  order  differential  of 
2z  sin  y  =  D.  There  results 

hy  -  -  Dm  •  Sz/2z2  cos  y  (A.  8) 

Furthermore,  if  v^  and^g  are  the  respective  velocity  vectors  of  the  planet  and 
the  spacecraft  and  if  6 r  is  the  difference  between  the  observed  and  the  reference 
occultation  times,  we  have 

Sz  =yp  Sr-  ( 5 r  +  vs  8r)  (A.9) 


=  -  S_r_  ~  yr  5  r 


where  v^  is  the  velocity  of  the  spacecraft  relative  to  the  planet. 


bining  Eqs.  (A.  7),  (A.  8) 

8t  - 


and  (A.  9)  we  have  finally 

1 


{p  -  tan  y  m)  •  vr 


(£-  tan  y  m)  •  8r 


Then  by  com- 
(A. 10 ) 


where p  is  a  unit  vector  perpendicular  to  m  and  lying  in  the  plane  determined  by 
the  lines-of- sight  to  the  planet  and  the  star. 


Star  Elevation  Measurement 

Consider  next  the  measurement  of  the  angle  between  the  lines-of- sight  to 
a  star  and  the  edge  of  a  planet  disc.  From  Fig.  A-2  we  have 

n  ‘  z  = z cos  (A  +  y)  ( A.  1 1 ) 
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where  A  is  the  angle  to  be  measured.  Again 

noting  that  6r  =  -  6z,  w'e  obtain 

— - — p  •  Sr  -  SA  +  Sy 
z  -  - 


taking  total  differentials  and 


(A.  12) 


-  SA  +  Dm  •  Sr/2z 


2 


cos  y 


=  S  A  +  tan  y  m  *  S_r  /  z 


or  finally 


S  A  -—[(>-  tan  y  m)  ‘  8_r 


(A. 13) 


Fig.  A-2.  Measurement  of  star  elevation  angle. 

Landmark  Measurement 

For  the  measurement  of  the  angle  between  a  landmark  on  a  planet  surface 
and  a  star,  let  p  be  a  unit  vector  perpendicular  to  the  line- of- sight  to  the  land¬ 
mark  and  in  the  plane  of  the  measurement.  Then  if  p  is  the  vector  position  of 
the  landmark  relative  to  the  center  of  the  planet,  we  have 

&  '  Sr, 


SA  =' 


U  +£\ 


(A. 14) 
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Radar  Range,  Azimuth,  and  Elevation  Measurements 

Assume  the  radar  site  to  be  the  origin  of  the  coordinate  system  although 
other  origins  could  equally  well  be  used.  Let  a  cartesian  coordinate  system  be 
chosen  such  that  the  z  axis  is  radially  out  from  the  center  of  the  Earth  through 
the  radar  site;  the  x  axis  is  positive  in  the  direction  from  which  radar  azimuths 
are  to  be  measured;  the  y  axis  completes  the  coordinate  system.  Then,  we  maj 

write  "  cos  0  cos  6 

cos  0  sin  6 
sin  0 


(A.  15! 


where  r,  9,  0  are,  respectively,  the  range,  azimuth,  and  elevation  of  the 
vehicle  from  the  radar  site.  Taking  differentials  separately  for  each  of  the 
three  variables  gives 

ii  cos  c  cos  a  ii 

(A.  16) 


Sr 


cos  0  cos  6 
cos  0  sin  6 
sin  0 


8  r 


Sr  =  r 


-  sin  0  cos  9 

-  sin  0  sin  6 

cos  0 


10 


(A.  17  ) 


8  r  -  i 


cos  0  sin  6 
cos  0  cos  6 
0 


(A.  18) 


S  9 


Then, 


by  expressing  each  of  these  relations  in  the  form  of  Eq. 
Sr  "  ||  cos  0  cos  6  cos  0  sin  9  sin  0  j|  S_r 

S/3  =-jJ-  ||  -  sin  0  cos  9  -  sin  0  sin  9  cos  0  ||  5_r 


S  6  = 


1 


r  cos  0 


sin  6  cos  6 


Sr 


(A.  1 ) ,  we  obtain 
(A.  19) 

(A. 20  ) 

(A. 21) 
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APPENDIX  B 


OPTIMUM  SELECTION  OF  NAVIGATION  MEASUREMENTS 


In  the  main  body  of  this  paper  a  method  of  processing  measurement  data 
in  an  optimum  linear  manner  has  been  developed.  The  purpose  of  this  appendix 
is  to  treat  the  associated  problem  of  selecting  those  measurements  which  are, 
in  some  sense,  most  effective.  For  example,  the  requirement  might  be  to 
select  the  measurement  to  be  made  at  time  t  in  order  to  get  the  maximum 
reduction  in  mean-squared  positional  or  velocity  uncertainty  at  time  t  .  Of 
perhaps  greater  significance  would  be  the  requirement  of  selecting  the  measure¬ 
ment  which  minimizes  the  uncertainty  in  any  linear  combination  of  position  and 
velocity  deviations.  Specifically,  one  might  select  the  measurement  which 
minimizes  the  uncertainty  in  the  required  velocity  correction.  As  a  further 
example,  one  might  wish  to  select  that  measurement  which,  if  followed  imme¬ 
diately  by  a  velocity  correction,  would  result  in  the  smallest  position  error 
at  the  target. 


Consider  first  the  simplest  case,  i.e.,  minimizing  the  mean-squared 
positional  uncertainty  at  time  t  .  From  Eq.  (4.  12)  the  mean-squared  positional 
uncertainty  is  expressible  as 


2 


e 


n 


(I)' 

n 


(B.  1) 


assuming  the  measurement  errors  to  be  uncorrelated.  In  the  absence  of  any 

2 

measurement  error  (  a  =  0  )  ,  the  problem  of  minimizing  either  mean-squared 

error  is  equivalent  to  finding  a  direction  for  the  h^  vector  which  maximizes  the 

ratio  of  two  quadratic  forms.  For  the  case  of  the  mean-squared  positional 

error,  the  geometrical  interpretation  is  clear.  Since  the  principal  directions 

of  E(  1 )  and  E^  are  the  same,  the  optimal  direction  for  h  coincides 

n  n  n  (jq  -n 

with  the  major  principal  direction  of  E' 
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The  problem  of  minimizing  the  mean-squared  velocity  uncertainty  at  time 
tn  by  proper  choice  of  the  h^  vector  is  not  as  easily  solved  or  interpreted. 
Again,  from  Eq.  (4.  12)  the  mean-squared  velocity  uncertainty  may  be  written 


'MED 


T  (2)'  p(3) '  , 
(4)  \  -n  En  En  On 


(B.  2) 


,T  pH)'  ,  .  2 

hn  tn  h„  +  an 


Denote  by  p  and  q  the  two  quadratic  forms 

P=hnTEf'Ei3)'hn  ,  q=hjEi,)'hn 


(B.  3) 


From  the  theory  of  quadratic  forms  there  exists  an  orthogonal  transformation 
which  will  reduce  q  to  a  diagonal  form.  Thus 


i?n 


(B.  4  ) 


gives 


T  nT  P(l)'  n  .  ,2  ,2  x  ,2 

q  -  d  U  En  L)  d  -  di  f  °2  +  ^3  “3 


(B.  5) 


(l)i 

where  H,,  /j.ot  are  the  characteristic  roots  of  the  matrix  E  and  the 

1  l  3  n 

columns  of  the  Q  matrix  are  the  associated  characteristic  unit  vectors.  Since 

(l)i 

E  ^  is  a  positive  definite  matrix,  the  characteristic  roots  are  positive  and  a 
further  transformation 


gives 


f  =  D  d 


q  =iT  i  =f?  +  f3 


(B.6 ) 


(B.  7) 


where  D  is  a  diagonal  matrix  whose  diagonal  elements  are  yf/7^,  yf  n0,  y/7i3- 

The  same  transformation  from  hn  to  f_  applied  to  the  quadratic  form  p 


produces 


p  =fT  D'1  QT  E<2)'  eJ,3)'  QD'1  f 


(B.  8) 


One  final  transformation  applied  to_f  will  reduce  Eq.  (B.  8)  to  a  diagonal  form 


results  in 


2  v  2  2 
p  -  A. i  mi  +  m 2  +  ^3  m3 


(B.  9  ) 
(B. 10 ) 


where  the  columns  of  the  S  matrix  are  the  characteristic  unit  vectors  of  the 
matrix  D  *  E  ^  1  E  ^  1  QD  *  and  A^,  >2,  A^,  the  corresponding 

characteristic  roots.  The  same  transformation  (B.9)  applied  to  (B.  7)  gives 


_  T  ,T  - 
q  -  _m  S  5  m 


:  m  ^  +  m2  +  m3 


(B.  11 


since  S  is  an  orthogonal  matrix. 


In  summary,  then,  the  transformation 


bn 


QD'1  S  m 


produces  for  the  ratio  of  the  two  quadratic  forms 

2  ,  ,  2 

A  n 

J> 

q 


Aj  m j  +  A>2  rn 2  ^3 


(B. 12) 
(B. 13) 


A  2 ) 


M 


rn o 


e(2) '  e(2) ’T 
n 


is  nonsingular,  the  product  E 


Furthermore,  if  the  matrix  E 

n  =  ‘  n 

is  positive  definite  and  it  would  then  follow  that  A 


(2)’  ^(3) 


n 


1' 


2’ 


are  all  real  and  positive. 

The  problem  of  maximizing  the  ratio  p/q  is  now  readily  solved.  Since 

no  measurement  error  is  assumed,  one  cannot  hope  to  determine  more  than 

the  direction  for  the  optimum  hn  or,  equivalently,  the  optimum  m.  Therefore, 

it  may  be  assumed  that  m.  is  a  unit  vector.  Let 

A^  =  max  f  A  j ,  A2 ,  A3) 

Then  the  optimum  value  of  m  is 

—  i  1  j  -  K 

(B. 15) 


(B. 14) 


"J 


irk 


The  same  technique  can  be  used  to  select  that  direction  for  h^  which 
minimizes  the  uncertainty  in  any  linear  combination  of  position  and  velocity 
deviations.  Specifically,  consider  the  selection  of  that  measurement  which 
minimizes  the  uncertainty  in  the  velocity  correction  which  would  be  required 
immediately  following  the  measurement. 


The  correlation  matrix  of  the  velocity  correction  uncertainty  is 

,T 


B  E  B 1 

n  n  n 


and  the  mean-squared  uncertainty  may  be  expressed  as 


*1  e;  bJ) 


h!  w  hn 

-  n  -n 


(B. 16) 


(B. 17) 


hT  E(1)'h  +  ~aA 

-n  un  -n  n 
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Here  W  is  a  symmetric  matrix  defined  by 


W  = 


P0)*  c(2)' 


bJ  Bn 


(B. 18) 


so  that  if 
definite. 


,(1) 


,(2) 


B  is  nonsingular,  the  matrix  W  will  be  positive 


n  n 

Under  any  circumstances,  if  the  identification 


B 


T 

n 


f  is  made,  then  the  exact  same  procedure  may  be  used  to  select  the  optimum  dir¬ 

ection  for  the  h^  vector  as  was  used  previously  to  minimize  the  mean-squared 
velocity  uncertainty. 


In  all  cases  of  practical  interest  the  determination  of  the  optimum  direction 
for  the  h^  vector  must  be  made  subject  to  certain  constraints.  For  example, 
one  might  wish  to  select  the  "best"  star  to  be  used  in  measuring  the  angle 
between  the  line  of  sight  to  the  center  of  a  planet  disc  and  the  line  of  sight  to 
the  star.  For  such  a  measurement  the  h  vector  is  required  to  be  perpendicular 
to  the  line  of  sight  to  the  planet.  If  z  is  the  position  vector  of  the  planet  from 
the  space  vehicle,  then  we  must  have 


hi  z  =  0 

-  n  -  n 


Applying  the  transformation  defined  in  Eq.  (B.  12  )  gives 

=  0 


T  CT  n- 1  nT 
m  S  D  Q  z„ 


(B.  19] 


(B. 20 ) 


T  -1  T 

Let  p  be  a  unit  vector  in  the  direction  of  S  D  Q  z  . 
—  —  n 


Then  the  problem 


of  selecting  the  optimum  direction  for  hn  or,  equivalently,  for  m  is  to  maximize 


2 

ml 


+  k> 


2 
m  ^ 


+ 


2 

m3 


subject  to  the  conditions  of  constraint 

p  =0  and  m  =1  (B.  21) 

In  terms  of  the  Lagrange  multipliers  p  and  a,  this  is  equivalent  to  the  problem 
of  obtaining  a  free  maximum  for 


E  *j  ■ 2p  E  pj  mj 

j=i  i=i 


2 

mj 


-  1 


I 

\ 
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Setting  the  partial  derivatives  with  respect  to  each  of  the  m.'s  equal  to  zero, 
we  have 


P  Pi 


I  Wj  -  a 


j  =1,  2,3 


(B. 22) 


where  p  and  o  are  to  be  determined  from  the  requirements  of  Eq.  (B-21). 


The  condition  that  m  be  orthogonal  to  p  leads  to  a  quadratic  equation 


for  cr. 


P]  (^2  +  ^3)  +  p2  ^1  +  *-3)  +  p3  + 


(B.  23) 


+  K>2  ^>2  P2  ^1  ^3  P3  ^2  ^ 

If  the  X  's  are  ordered  X^  <  ^2  <  ^3'  t*ien  two  roo'ts  CT1  and  CT2 
be  such  that  X1<a1<X2<  a2<  Xg_  Qne  of  these  roots  provides  the 

maximum  while  the  other  gives  the  minimum.  The  other  Lagrange  multiplier 
p  is  determined  so  that  m  will  be  a  unit  vector. 


With  the  optimum  vector  m  selected,  the  corresponding  value  for  h^  is 
found  from  Eq.  (B.  12)  . 
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